load datos3b.mat
	
for i =1:length(xd)
	if xd(i) != 0
		x(i) = 1./xd(i);
	else
		x(i) = xd(i);
	end		
	if yd(i) != 0
		y(i) = 1./yd(i);
	else
		y(i) = yd(i);
	end		
end

P = polyfit(x,y,1)

A = P(1) % la ordenanda al origen
B = P(2) % la pendiente

		
plot(xd,yd,'r',xd,xd./(A*xd + B),'b')

